Suppressing noise-induced intensity pulsations in semiconductor lasers by means of 

time-delayed feedback 
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We investigate the possibility to suppress noise-induced intensity pulsations (relaxation oscilla- 
tions) in semiconductor lasers by means of a time-delayed feedback control scheme. This idea is 
first studied in a generic normal form model, where we derive an analytic expression for the mean 
amplitude of the oscillations and demonstrate that it can be strongly modulated by varying the 
delay time. We then investigate the control scheme analytically and numerically in a laser model of 
Lang-Kobayashi type and show that relaxation oscillations excited by noise can be very efficiently 
suppressed via feedback from a Fabry-Perot resonator. 



I. INTRODUCTION 

In many dynamical systems noise plays an important 
role and influences the system's properties and the dy- 
namic behavior in a crucial way. Control of the noise- 
mediated dynamic properties is a central issue in nonlin- 
ear science 

An often encountered effect of noise is the excitation 
of irregular stochastic oscillations under conditions where 
the deterministic system would rest in a stable steady 
state, e.g., a stable focus. The random fluctuations then 
push the system out of the steady state. These noise- 
induced oscillations are a widespreadphenomenon and 
appear, for instance, in lasers [2, 3, 4 S 01, chemical re- 
action systems |7], semiconductor devices neurons 
and many other systems. 

In practical applications the need arises to control the 
oscillations, for instance, by increasing their coherence 
and thus the regularity of the oscillations. In recent years 
different methods to control stochastic systems have been 
developed, and applied to noise-induced oscillations in 
a pendulum with a randomly vibrating suspension axis 
and external periodic forcing [11], stochastic resonance 
[H, [Ts*], noise-induced dynamics in bistable delayed sys- 
tems [14, 15], and self-oscillations in the presence of noise 
In the context of coherence resonance [13, fl8| . 
time-delayed feedback in the form originally suggested 
by Pyragas to stabilize unstable states in deterministic 
systems [l|, [l§| has been demonstrated to be a power- 
ful tool to control purely noise-induced oscillations (2Q| . 
This method couples the difference of the actual state 
X(t) G R" and of a delayed state of the system X[t — r) 
back into the system. 



dt 



X{t) = f{X{tlt) - K[X{t)-X{t-r)]. 



The time delay r and the (matrix valued) control ampli- 
tude K are the control parameters which can be tuned. 
While previous studies have shown that the delayed feed- 
back method can control the main frequency and the cor- 
relation time tcor and thus the regularity of noise-induced 
oscillations in simple systems [20, 21, 22, 23^4, 25, 26] 
as well as in spatially extended systems [23, , and 

deteriorate or enhance stochastic synchronization of cou- 



pled systems [30], in this paper we focus on the suppres- 
sion of stochastic oscillations. We analyse the mean am- 
plitude (or, more generally, the covariance) of the oscil- 
lations and show that time-delayed feedback control can 
decrease the mean oscillation amplitude for appropriately 
chosen delay time, and thus suppress the oscillations. 

The paper is organized as follows. In section [TTl we 
study a generic model consisting of a damped harmonic 
oscillator driven by white noise and investigate the influ- 
ence of delayed feedback. In this generic system we de- 
rive an analytic expression for the mean square oscillation 
amplitude in dependence on the feedback and show how 
the oscillations can be suppressed. In section [TTll we con- 
sider a semiconductor laser, a practically relevant exam- 
ple, and show how optical feedback from a Fabry-Perot 
resonator, which realizes the delayed-feedback scheme, 
can suppress noise-induced relaxation oscillations in the 
laser. 



II. GENERIC MODEL 

We consider a damped harmonic oscillator (whose flxed 
point is a stable focus) subject to noise {^) and feedback 
control 



z(t) = {X-iLoo)z{t) + Dat) 
-K[zit)-z{t-r)l 



(z e C) (1) 



where A < and coq are the damping rate and the natural 
frequency of the oscillator, respectively, D is the noise 
amplitude, K is the (scalar) feedback strength and r is 
the delay time of the control term. We consider Gaussian 
white noise 

m = + '-6(0, (e- e R) 
(e-) = 0, 
(e-(t)o(f')) = s,j5it-t'). 

In our particular system the delay term in eq.([l]) does not 
induce any local bifurcations in the deterministic system. 
Thus, the flxed point is stable for all r and K. 

A similar normal form (without noise) was previously 
used to study the stabilization of unstable determinis- 
tic flxed points by time-delayed feedback [3l|, which is 



possible in the same way as stabilization of unstable de- 
terministic periodic orbits [J, [H, [13, [Hi . 

The power spectral density of z has been calculated in 
(23I and is given by (see Fig. [2]) 



27r [A - - cos(cjr)]2 [uj - ujq ^ K s\n{u;r)]^ ' 



Figs. [T] and [2] display the dependence of the power spec- 
tral density on the delay time. With increasing delay new 
peaks appear in the spectrum. These are related to new 
modes generated by the delay as will be shown later. 




Figure 1: (color online) Power spectral density 5 as a function 
of the frequency co and the delay time r (logarithmic color 
scale). Parameters: A — —0.01, cjo = 1 (To = 27r), D = 1, 
K = 0.2 
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Figure 2: (color online) Power spectral density 5{ou) for dif- 
ferent delay times r. Parameters: A = —0.01, cuq = 1, D = 1, 
K = 0.2 

Before proceeding, we transform eq.([l]) into a rotating 



frame z(t) = u{t) e''^^^ 

b{t) = {X - K) u{t) ^ K e'""'^ u{t - t) ^ e'"^'' D ^{t) 
= au{t)^bu{t-r)^Di{t), (2) 

where a G M, 6 G C, and = e''^o^^(t) is a noise term 
with the same properties as ^(t). The purpose of the 
transformation is to make the parameter a real, which 
will be necessary later. 

In [35I Kiichler and Mensch analyzed equation ([2j) for 
real variables. We will follow their approach and adapt 
it to complex variables. Similar results for the Van der 
Pol oscillator have been obtained independently in (36| . 
A different two-dimensional system with noise and delay 
has been recently studied in [37]. 

We will calculate the autocorrelation function 



G{t) = {u{s + t)u{s)) 

in an interval t e [0,r], where the overbar denotes com- 
plex conjugate. In particular, this gives the mean square 



amplitude (r 



= G(0) of the oscilla- 



tions. With the Green's function uo{t) solving 

uo{t) - a uo{t) - buo{t - r) = S{t), 

with iyo(t) = for t < 0, we can formally find a solution 
of equation ([T]) 



u{t) = / dtiUo{t-ti)Di{ti). 



(3) 



Using (|3]) we obtain 
G(t) = {u{t + 1)^)) 

t+t t 



J dti J dt2 uoCt + t - ti) uoCt - t2) 



-00 —00 



2D^ J ds uo{s ^ t) uo{s) 


= 2D^ C{t). 

The Green's function uq can be calculated [H, [3§| by 
iteratively integrating eq. ([2]) on intervals [kr, (/cH- l)r) 



k=0 



From the definition of C and uq it follows that C satisfies 
the following equations 



C(t) = C(-t) (4) 
C{t) = aC{t) + bC{t-T) (t>0) (5) 
C{t) = aC{t) + bC{T-t) (t>0). (6) 



Using these three equations, we can find an ordinary dif- 
ferential equation for C, using eqs.®,©, 

c/2 



C{t) = aC'{t)- bC{T-t) 
= a[aC{t) + bC{t-T)] 



-b[aC{T-t) + bC{-t)] 
a^C{t) + abC{t-T) 



-abCit-T)-\b\^Ci-t) 
= {a'-\bf)C{t). 

Here it was necessary to have a real a, in order for the 
delay terms to cancel. Thus C is of the form 



C{t) = Ae^'^Be 



-At 



with 



The complex coefficients A and B can be found from the 
equations 



and 



C(0) = C(0) G M, 
C(0) = aC(0) + foC(7) 



oo 

j c/s^[(7o(s)tVo(s)] 


oo 

j ds [uo{s) uo{s) + uo{s) uo{s)] 



(7) 

(8) 



(9) 



= a C(0) + b C(r) + a C(0) + 6 C(r). 

Solving equations (O,® and (|9]) for A and 8 gives the 
mean square oscillation amplitude 



C(0) = (r^) = Re(/\) + Re(8) 
1 



/<2 + - /<2cosh(2Ar) 



4A K cosh(Ar) [Acos(cJor) + K sinh(Ar)] + a [A + K cos{u;or) sinh(Ar)] " 

I 



(10) 




100 



Figure 3: (color online) Mean square amplitude (r^) of noise- 
induced oscillations (solid line: analytics, dots: numerics). 
Parameters: A = -0.01, D = 1, K = 0.2, cjq = 1 (To = 27r) 
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Figure 4: (color online) Mean square amplitude (r^) of oscilla- 
tions as a function of the delay time r for different K (analytic 
solution). Parameters: A = —0.01, D = 1, cjo = 1 (7~o = 27r) 



This is now an analytic result which allows to analyze 
the effect of the control term. Figure [3] displays analytic 
and numeric results for the oscillation amplitude. The 
dependence on the control force K is shown in Fig. [H 

The oscillation amplitude can thus be strongly modu- 
lated by varying r. We can obtain the envelopes of the 
modulation by setting the terms cos(a;oT) to their maxi- 



mum and minimum values ±1 in eq. (p!Ql 



/<sinh(Ar) =F A 



2A /<cosh(Ar)±a' 

Figure [Ifa) displays (r^) and the envelopes versus r. 
The mean square oscillation amplitude is modulated as 
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Figure 5: (color online) (a) Mean square amplitude (r^) of 
oscillations and (b) real part of the eigenvalue spectrum of 
the fixed point as a function of the delay time r. The dashed 
lines in (a) mark the envelope C±. Parameters: A = —0.01, 
D = 1, K = 0.2, Lvo = 1 (To = 27t) 



eigenvalues corresponds to the damping rate of the re- 
spective mode, the oscillation amplitude excited by noise 
is large if a mode is weakly damped, and small if all 
modes have rather large (negative) damping rates. 

Because eq. ([T]) is linear and ^ is Gaussian noise, the 
probability distribution pjx y), where z = x + /y, is also 
a Gaussian distribution [33|. The rotational invariance 
{z' = z e'^) of eq. ([T]) implies, that p(x,y) is invariant 
under rotations, too. These two arguments lead to the 
probability distribution 



P(x,y) = - 




exp 



with 



Figure [6] shows the marginal distribution 



dyp(x,y) 



a function of r with a period Tq = 27r/a;o. The maxima 
and minima occur at 

r+ = A? To and r_ = — - — Tq 

respectively. The smallest oscillation amplitude is 
reached at 



I opt 



To/2. 



To understand the behavior of the mean square oscil- 
lation amplitude as a function of the delay time r, one 
has to look at the eigenvalue spectrum of the fixed point 
z = of eq. ^ (without fiuctuations). The ansatz 
z(t) (X e^^ in eq. ([T]) gives rise to a transcendental equa- 
tion for the eigenvalues /i: 



fi = {X- iuo- K)^ Ke-^^. 



(11) 



This equation can be solved using the Lambert function. 
The Lambert function W is defined j[39|] as the inverse 
l/l/(z) of the equation 



We 



w 



{ze 



(12) 



Since eq. ([T2|) has infinitely many solutions, the Lambert 
function W has infinitely many branches l/l/n(z) indexed 
by n. Using the Lambert function W the solutions of ([TT]) 
are given by 



r^n = WnirKe 



-{X—iooo — K)i 



] + (A - icoo - K)t. 



FiglH^b) shows the real part of the spectrum versus r. 
As r increases, different eigenvalue branches originating 
from — oo approach the zero axis (albeit remaining < 0) 
and then bend away again. Since the real part of the 



for r = (dashed) and r = Topt = To/2 (solid). The case 
r = corresponds to no control, because the feedback 
term in ([T]) vanishes. The case r = Topt realizes the opti- 
mal delay time, where the oscillations are most strongly 
suppressed and the distribution p(x,y) is narrowest. 
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Figure 6: (color online) Marginal probability distribution p(x) 
without (r = 0) and with optimal control (r = Tq/2) Param- 
eters: A -0.01, K 0.2, D l,cjo = 1 (To = 27r) 



III. LASER MODEL 

In this section we investigate the effects of feedback 
and noise in a semiconductor laser. A laser with feedback 
from a conventional mirror can be described by the Lang- 
Kobayashi equations lioll. Other types of feedback have 
also been investigated (4lL[42|. One particular feedback 
realizes the delayed feedback control with an all-optical 
scheme [13, (H]. The feedback is here generated by a 



E complex field amplitude 

n carrier density 

a linewidth enhancement factor 

K feedback strength 

r roundtrip time in the Fabry-Perot 

p excess pump injection current 

T timescale parameter 

Fe noise term describing spontaneous emission 

P spontaneous emission factor 

no threshold carrier density 

(f, il^ phases depending on the mirror positions 

Table I: Nondimensionalized laser variables and parameters 



Fabry-Perot resonator. A schematic view of this setup 
is shown in Fig. [3 A fraction of the emitted laser light 
is coupled into a resonator. The resonator then feeds an 
interference signal of the actual electric field E{t) and the 
delayed (by the round trip time) electric field E{t — r) 
back into the laser. 
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Figure 7: (color online) Setup of a laser coupled to a Fabry- 
Perot resonator realizing the time-delayed feedback control 

A modified set of non-dimensionalized [i^ Lang- 
Kobayashi equations |;42i] describes this setup 



dt 



-{l + ia)nE (13) 
-e''^K[E{t)-e''^E{t-T)] + FE{t), 



Tj^n = p-n-{l + n)\Ef, 

where the variables and parameters are defined in Table 

n 

The phases (p and i/j depend on the sub-wavelength 
positioning of the mirrors. By precise tuning cp = 27rn 
and tjj = 27r/T7 one can realize the usual Pyragas feedback 
control 

-K[Eit)-E{t-T)]. 

We consider small feedback strength K, so that the laser 
is not destabilized and no delay-induced bifurcations oc- 
cur. A sufficient condition 1431 is that 



K <Kc 



rVTTo^ 



The noise term Fe in (p!3|) arises from spontaneous emis- 
sion, and we assume the noise to be white and Gaussian 



with the spontaneous emission rate 

Rsp = f3{n^no). 

where f3 is the spontaneous emission factor, and a7o is 
the threshold carrier density. Without noise the laser 
operates in a steady state (continuous wave cw emission) . 
To find these steady state values, we transform eqs. (p!3|) 
into equations for intensity / and phase (j) hy E = ^/l e'^ 
(see Appendix A): 

/ = A7/-2/<[/-^V^cos((/), -(/))] + /?sp + F/(t), 



dt 

±^ = ]^an^K^ sin((/>. -</>) + F^t), (14) 
T—n = p- A7- (1 + A?)/, 
where Ir = l{t — r), (j)r = (j){t — r), and 

(F,)=0, (F^)=a 

(F,(t)F^(0) =0, 

(F,(t)F,(0) = 2RsplS{t-t') 

Setting Tt' = Tt^ = ft^ = ^o"^^' /< = and 
replacing the noise terms by their mean values, gives 
a set of equations for the mean steady state solutions 
A?* and (j) = uo^t without feedback (the solitary laser 
mode) . Our aim is now to analyze the stability (damping 
rate) of the steady state. A high stability of the steady 
state, corresponding to a large damping rate, will give rise 
to small-amplitude noise-induced relaxation oscillations 
whereas a less stable steady state gives rise to stronger 
relaxation oscillations. Linearizing eqs. (p!4|) around the 
steady state X{t) = X* + (5X(t), with X{t) = {I, (j), n) 
gives 



dt 

with 



X{t) = UX{t)-V [X{t) -X{t-T)] + F{t), (15) 



U = 



I* +13 



- f; 


^l+n,) -^(1+/*) 



(Fe) = 0, {FE{t) FE{t')) = RspS{t - t'). 



V = d\ag{K, K, 0) 
where diag{...) denotes a 3 x 3 diagonal matrix, and 

F = {Fi, F^, 0). 
The Fourier transform of eq. ([T5|l gives 

X{oj) = [iu}-U+V{l- e"'"^)]-^ 

^ V ' 

=M 



The Fourier transformed covariance matrix of the noise 
is 

{F{u,)F{uy) = i-diag(2/?,p/., 0)5{u>-u>'), 

with the adjoint f. The matrix- valued power spectral 
density S(cj) can then be defined through 

and is thus given by 

S{uj) = 6\ag{Ssi{uj), Ss(i,{uj), Ssn{^)) 
= ^Md\ag{2RspU, 0)Mt. 

The frequency power spectrum is related to the phase 
power spectrum Ss(f){u;) by [Hi 




Figure 8: (color online) Analytical (left) and numerical (right) 
results for the power spectral density Ssi{io) of the intensity 
for different values of the delay time r. 

Parameters: p = 1, T = 1000, a = 2, (3 = 10"^ no = 10, K = 
0.002. (A typical unit of time is the photon lifetime Tp = 
10~^^s, corresponding to a frequency of 100 GHz.) 



analytics 





0.01 
0.008 
200 0.006 
150 0.004 

00 ^-^^^ 



Figure 9: (color online) Analytical (left) and numerical (right) 
results for the power spectral density Ss(j){oo) of the frequency 
for different values of the delay time r. 

Parameters: p = 1, T = 1000, a = 2, f3 = 10"^ no = 10, K = 
0.002 



the nonlinear simulations. The main peak decreases with 
increasing r and reaches a minimum at 



Tro 



2lT 



I opt 



2Q 



100. 



RO 



With further increasing r the peak height increases again 
until it reaches approximately its original maximum at 
T ^ Tro. a small peak in the power spectra indicates 
that the relaxation oscillations are strongly damped. 
This means that the fiuctuations around the steady state 
values U and a?* are small. Figure [10] displays exemplary 
time series of the intensity with and without feedback. 
The time series with feedback show much less pronounced 
stochastic fiuctuations. 
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Figure 10: (color online) Intensity time series with (top panel) 
and without (bottom panel) control 



Parameters: p = 1, T = 1000, a = 2, /3 = 10" 
100 ^ ro/2 



no = 10, r = 



Next, we study the variance of the intensity distribu- 
tion as a measure for the oscillation amplitude 



(/-(/)r 



Figures [8] and [9] display the intensity and the frequency 
power spectra, respectively, for different values of the 
delay time r, obtained analytically from the linearized 
equations (left) and from simulations of the full nonlin- 
ear equations (right). All spectra have a main peak at the 
relaxation oscillation frequency Q.ro ^ 0.03. The higher 
harmonics can also be seen in the spectra obtained from 



This measure corresponds to the quantity (r^) which we 
have considered in Section IL Figure [TT] displays the vari- 
ance as a function of the delay time. The variance is min- 
imum at r ~ Tro/ 2^ thus for this value of r the intensity 
is most steady and relaxation oscillations excited by noise 
have a small amplitude. This resembles the behavior of 
the generic model (see Fig. [5](a)). 
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Figure 11: (color online) Variance of the intensity / vs. the 
delay time. Parameters: p = 1, T = 1000, a = 2, /3 = 
10"^ no = 10, K = 0.002 

Figure [12] displays the intensity distribution of the laser 
without (dashed) and with (solid) optimal control (com- 
pare Fig. [6]). The time-delayed feedback control leads to 
a narrower distribution and less fluctuations. 
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Figure 12: (color online) Probability distribution of the inten- 
sity / with and without the resonator (simulations). Param- 
eters: p = 1, r = 1000, a = 2, /3 = 10-^ no = 10, K = 0.002 



IV. CONCLUSION 

In this paper we have shown that time-delayed feed- 
back can suppress noise-induced oscillations. 

In the flrst part we investigated a generic normal form 
model consisting of a stable focus subject to noise and 
control. We found an analytic expression for the mean 
square amplitude of the oscillations. This quantity is 
modulated with a period of Tq = Itt/ujo in dependence 
on r. For r = To/2 the oscillations have the smallest 
amplitude. 

In the second part we considered a semiconductor laser 
coupled to a Fabry-Perot resonator. In the laser spon- 



taneous emission noise excites stochastic relaxation os- 
cillations. By tuning the cavity round trip time to half 
the relaxation oscillation period Topt ~ Tro/2 the oscil- 
lations can be suppressed to a remarkable degree. This 
is demonstrated in the power spectra of the intensity and 
the frequency, where the relaxation oscillation peak has 
a minimum height at Topt- The variance of the inten- 
sity distribution A I shows a minimum at Topt, thus the 
intensity distribution is narrowest at this value of r. 
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Appendix A: ITO TRANSFORMATION 

Ito's formula describes how a stochastic differential 
equation (SDE) is transformed to new coordinates. Con- 
sider the stochastic differential equation for x(t) 

dx{t) = a[x(t), t] dt + b[x{t), t] dW{t). 

Ito's formula specifles the transformation to a new vari- 
able y = f{x). The SDE for y is given by 

dy = df[x{t)] 

= a[x(t), t] f'[x{t)] dt + b[x{t), t] f'[x{t). t] dW{t) 

+ ^b[x{t),tff"[x{t)]dW'. 

We will apply Ito's formula to rewrite the laser equa- 
tions for the complex electric fleld E in terms of the am- 
plitude A and the phase (j) 

The equation for E without feedback is given by 

±E=\{l + ia)nE + FE{t) 
or written as a stochastic differential equation 

c/E = ^(1 + ia) nEdt^ dW{t), 

with the complex Wiener process dW = dWx + i dWy. 
We deflne the new coordinates 

/i + /(^ = log/\ + loge'^ = log[Ex + /^y]. 

Using Ito's formula with 

a[E] = ^{l+ia)nE. b[E] = ^J^ and, f[E] = log E. 



we find 



with 



d{^^+i4>) = lil+ia) n 1 c/t-i^^ dW'+^ ^ ^ dW. 

(Al) 

For a complex Wiener process dW one can easily see that 
dW^ = [dW^ + idWyf = d\Nl + 2/ d\N^d\Ny - dW^ = 0. 

Here we used dW^ = dW^ = dt and dW^dWy = HJ. 
Thus, equation (|Aip simplifies to 



d{fi-\- = ^ '^^ ^ exp{—fi — dt 



+ 1 



exp(— /i — dW . 



Splitting this equation into real and imaginary part and 
transforming with Ito's formula back to A = exp/i, we 
obtain 

dA = {^n A^^)dt^]J^ {cos (I) dW^^sm (I) dWy) 

1 1 rR 

dcj) = -A?acft+-A/^(-sin0c/l/14 + cos^c/l/l/y). 

Because the rotation is an orthogonal transformation, 
one can understand the increments as new independent 
Wiener processes 

dWA = cos^c/l/l/x + sin0c/l/l/y 
dW^ = -sin(/)c/l/l/x + cos(/)c/l/l/y. 

We have derived the laser equations in polar coordi- 
nates. To include the delay terms does not change the 
derivation and we will just state the result here 

^^A = ^nA - K[A - Ar cos((/., _</,)] + ^ + F^it) 



— = -an ■ 
dt 2 



{FA{t)F4t')) 



S{t-t') 



{F4t)F4t')) = ^sit-n 

To obtain the equations for intensity / = f{A) = A^ 
instead of the amplitude Ito's formula has to be applied 
again (of course this could be done in one step from the 
initial equations). The amplitude equation is given by 



^nA - K[A - Ar cos{(t)r - (^)] + !> c/t 




dA 



dW[t). 

For a real stochastic process dW holds dW^ = dt. Using 
Ito's formula with dW^ = dt, r{A) = 2A and, r{A) = 2, 
we find 

dl = {l-nA- K[A-ArCos{(t)r-(t))]^^)2Adt 
+ i/?,p2c/t+y^2/\dl/l/(t) 



and thus 



dt 



with 



/ = n\ -K\\ -\r\\JTrCO'-.{4,r-4>)\ 

+ Rsp + Fi{t), 
{Fi{t)Fi{t'))=2RsplS{t-t'). 
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